Quantum correlations in the thermodynamic limit: the XY-model 
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Abstract 

We investigate thermal properties of quantum correlations in the thermodynamic limit with 
reference to the XY-model 
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I. INTRODUCTION 



Since the formalization by Werner [T] of the modern concept of quantum entanglement it 
has become clear that there exist entangled states that comply with all Bell inequalities 
(BI). This entails that non-locality, associated to BI- violation, constitutes a non-classicality 
manifestation exhibited only by just a subset of the full set of states endowed with quantum 
correlations. Later work by Zurek and Ollivier j2] established that not even entanglement 
captures all aspects of quantum correlations. These authors introduced an information- 
theoretical measure, quantum discord, that corresponds to a new facet of the "quantumness" 
that arises even for non-entangled states. Indeed, it turned out that the vast majority of 
quantum states exhibit a finite amount of quantum discord. 

The tripod non-locality-entanglement-quantum discord is of obvious interest and possesses 
technological implications. The crucial role played by quantum entanglement in quantum 
information technologies is well known [3]. In some cases, however, entangled states are 
useful to solve a problem if and only if they violate a Bell inequality jl]. Moreover, there 
are important instances of non-classical information tasks that are based directly upon non- 
locality, with no explicit reference to the quantum mechanical formalism or to the associated 
concept of entanglement [5]. Last, but certainly not least, recent research indicates that 
quantum discord is also a valuable resource for the implementation of non-classical informa- 
tion processing protocols [6HTU] . On the light of these developments, it becomes imperative 
to conduct a systematic exploration of the connections between the tripod members. 
We investigate here the relation between quantum discord and entanglement and in an 
infinite system, namely, the XY model in the thermodynamic limit [TT]. This model, like 
the celebrated Ising and Heisenberg models, is one of the paradigmatic systems in statistical 
mechanics. The Hamiltonian of the anisotropic one-dimensional spin-^ XY model in a 
transverse magnetic field h (N particles) reads 

N N 

H = £[(1 + ^SiSt 1 + (1 - 1)SISI +1 ] -hJ^Si, (1) 
i=i j=l 

where a{ = 2S 3 U (u = x,y,z) are the Pauli spin-^ operators on site j, 7 G [0,1] and 
a i +N = °u- The model Q for iV = 00 is completely solved by applying a Jordan- Wigner 
transformation [TT| H5], which maps the Pauli (spin 1/2) algebra into canonical (spinless) 
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fermions. The system (except for the isotropic case 7 = 0) undergoes a paramagnetic-to- 
ferromagnetic quantum phase transition (QPT) [TJ1 US] driven by the parameter h at h c — 1 
and T = 0. It is well known that near factorization a characteristic length scale naturally 
emerges in the system, which is specifically related with the entanglement properties and 
diverges at the critical point of the fully isotropic model [16J. 

It is thus our intention in this communication to study the interplay of entanglement and 
quantum discord for the XY model in the thermodynamic limit. To such an effect we will 
consider in Section II the correlations existing between a pair of qubits located at two given 
sites. For comparison purposes we shall also discuss in Section III the correlations between 
pairs of qubits in the finite Heisemberg model. Some conclusions are drawn in Section IV. 

A. Quantum discord 

Quantum discord [21 16] constitutes a quantitative measure of the "non-classicality" of bi- 
partite correlations as given by the discrepancy between the quantum counterparts of two 
classically equivalent expressions for the mutual information. More precisely, quantum dis- 
cord is defined as the difference between two ways of expressing (quantum mechanically) 
such an important entropic quantifier. Let p represent a state of a bipartite quantum sys- 
tem consisting of two subsystems A and B. If S(p) stands for the von Neumann entropy of 
matrix p and pa amd pB are the reduced ("marginal") density matrices describing the two 
subsystems, the quantum mutual information (QMI) M q reads [2] 

M q (p) = S( PA ) + S( PB )-S(p). (2) 

This quantity is to be compared to another quantity M q (p), expressed using conditional 
entropies, that classically coincides with the mutual information. To define M q (p) we need 
first to consider the notion of conditional entropy. If a complete projective measurement IP? 
is performed on B and (i) pi stands for Ttab nf p and (ii) Pa\\u b f° r Pf pllf /pi], then the 
conditional entropy becomes 

S(A\{nf}) = Y,P>S(p Am ?), (3) 

i 

and M q (p) adopts the appearance 
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quantum discord A, that quantifies non-classical correlations in a quantum system, including 
those not captured by entanglement. One notes then that only states with zero A may exhibit 
strictly classical correlations. Among many valuable discord-related works we just mention 
two at this point that are intimately related to the present one, e.g., those of Zambrini et 
al. pi] and Batle et al. [E]. 

II. TWO QUBITS IN THE INFINITE XY MODEL 

The general two-site density matrix is expressed as 



R = j — i is the distance between spins, {u,v} denote any index of {<To, cr x , o~ y , a z }, and 
TmP = (o~ l u (g) a 3 v ). Due to symmetry considerations, only {Txx ,Tyy,Tzz ,Txy} do not 
vanish. Barouch et al [13] have provided exact expressions for two-point quantum corre- 
lations, together with details of the dynamics associated with an external field h(t). For 
the purposes of this paper, we will consider only systems which at time t — are in 
thermal equilibrium at temperature T. We have then the canonical ensemble expression 
pit = 0) = exp [— (3 H], where /3 = 1/kT and k is the Boltzmann constant. Following [1 3j , 



one obtains T^J = Tyy = Gi,tH& = — G\G-i — SiS-\, and T^y = Si, where 






(h - h f )(hf - cos0) cos(2A(h f )t)] 




x [{7 2 sin 2 + (h — cos 4>)(hf — cos (f>) } (cos (j> 




(6) 
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s « = — y„ #smw > sm ^ a(ma(m ■ (7) 



7T 



with A(/i) = [7 2 sin 2 + (h — cos 0) 2 ] 1 / 2 . G^ is the two-point correlator appearing in the 
pertinent Wick-calculations and M z = \Gq. The two-spin correlation functions are given by 



y y 



Gi Go 
G2 Gi 



Gr Gr-i 
K< +R ) = 4(^) 2 - G*G_ 



G-_r 
■■ G_! 

G-R+2 
G-R+3 

Gi 



(8) 



(9) 



(10) 



where R = j — i (distance between spins). In the case where more than two particles 
are considered, the previous correlators no longer possess their previous Toeplitz matrix 
structure [19] . 

It will prove convenient to cast the two qubit states (5) in two forms. States p\j are 
written in the computational basis {|00), |01), |10), |11)} as 



1 

4 



I 1 + 4M Z + T ZZ 



\ T xx Tyy ~\- i"2T X y 
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l-AM z + T zz ) 



(11) 



These very states p \a acquire instead the following form in the Bell basis {|$ + ), |$ ), |^ + ), |^ }} 



/ p xx ip[ 2 ip{ 3 pf 4 \ 

-ip[ 2 P22 P23 ipL 

-iPis Pis P33 ip T M 

\ PU -*P34 P44 / 



(12) 



This special form is such that one can use it to analytically compute the maximal violation of 
a Bell inequality, a measure for nonlocality [Tj5]. In turn states in (11) are of such special 
aspect that the quantum discord Qd turns out to be to be analytically given (see Ref. [2U]). 
Nevertheless, the concomitant Qd can be easily obtained, in different fashion, as follows. 
The most general parameterization of the local measurement that can be implemented on 
one qubit (let us call it B) is of the form {IL% = I A ® |0')(0'|,n^ = I A ® More 
specifically we have 



|0') <- cosa|0) + e i/3 ' sina|l) 

|T) 4- e-^sinalO) -cosa|l), (13) 

which is obviously a unitary transformation -rotation in the Bloch sphere defined by angles 
(a, /?')- for the B basis {|0), |1)} in the range a G [0, tt] and (3' 6 [0, 27r). After some 
cumbersome calculations, it turns out that the expression for a minimum discord A of 
the Introduction exhibits a positive and nonsingular Hessian, convex for the relevant range 
of values of (a, /?'). Our expression possesses thus a unique global minimum, that occurs 
when the concomitant partial derivatives vanish. This happens whenever we have (sin a = 
^, sin/3' = 0). 

The present results correspond to pairwise entanglement and quantum discord for the infinite 
XY model at any temperature, including zero-one. This implies that one does not really 
need to "solve" the model in the sense of sufficiently augmenting the number of spins in the 
chain for the results to be thermally relevant. T here is an actual, thermometer-measurable 
temperature, since we are tackling a "real" thermodynamic system. This is to be confronted 
to the vast XY-literature associated to finite spin-numbers, where T is not, strictly speaking, 
well-defined in the thermodynamics sense. 

A comparison between the discord Qd and the entanglement of formation E at T = is 
displayed in Fig. 1 (from now on we shall take the Boltzmann constant k — 1). Qd and 
E are depicted versus the external magnetic field h (anisotropy 7 = |) for the nearest 
neighbor configuration R = 1. Remarkably enough, the Qd measure exhibits a maximum 
in the vicinity of the factorizing field hf — y/l — 7 2 . Both Qd and E seem to decay in the 
same fashion. The classical correlations (CC) for the same configuration are depicted in the 
inset of Fig. 1. Notice that all quantities here considered, i.e., Qd, E, or CC, are ultimately 
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described in terms of several Grs for all configurations, so that they all diverge at the QPT 
(for h — 1) in the same way. 

As an illustration consider the magnetization given by M z (h) = \Gq = ^ lo d^Yf 2 sm2 4>+ 
(h - cos0) 2 ]V 2 . For 7 = 1 we have M z (h) = f h (^E[^]) = + 
^^-^(f+f)]) where K(E) is the complete elliptic integral of the first(second) kind. Since 
J^M 2 diverges in logarithmic way at h = 1, as also do the divergence of K and the first 
derivatives of E, Qd, and CC. In other words, they all signal the presence of a h = 1-QPT at 
zero temperature (except for the isotropic case 7 = 0). In fact, the possibility of detecting 
a QPT at finite T by using Qd has been recently considered by Werlang et al. [21]. They 
perfom an interesting analysis of the role of the temperature and Qd in several quantum 
systems. We remember that a different concept such as nonlocality -as measured by the 
maximum violation of the well known Clauser-Horne-Shimony-Holt Bell inequality- was also 
considered as a QPT in [Tj5] (also in the context of the XY-model). In the present work 
we do not focus attention on this particular issue of QPTs, but study instead the compari- 
son between entanglement and quantum discord for finite and infinite systems at non-zero 
temperature. 

Fig. 2 depicts the the same quantities as Fig. 1 for several configurations. As we increase 
the relative distance from R = 1 to 2, 3, and 00, the corresponding Qd's diminishes and also 
decays in faster and faster fashion with h. Notice that while entanglement (not shown here) 
globally diminishes, Qd only tends to vanish for h > 1 and R = 00. The inset here depicts 
the CC for the same configurations. They decreasing in the same fashion. While E tends to 
zero, both Qd and CC remain nonzero, regardless of the distance between spins along the 
infinite chain. 

As soon as we introduce a non-zero temperature things drastically change. In Fig. 3 we 
display several quantities at different temperatures (T = 0.01,0.1,0.3,0.5,1): R = 1 and 
7=|- Fig- 3(a) depicts the entanglement of formation E for states p\j (5) versus the 
magnetic field h as we increase the temperature. T lowers and broadens the region of null 
entanglement from a point at the factorizing field hf (T = 0) to finite intervals centered 
at hf. Eventually, E becomes finite at higher values of h. This temperature-generated 
entanglement is depicted quantitatively in Fig. 3(b), where the region of zero entanglement 
extends from a point at zero T to a finite-sized region as T grows. The aforementioned 
region ceases to be finite beyond a critical temperature that depends on the particular 



i?'s and 7's involved therein. We discern some resemblance with a phase diagram: within 
the area encompassed by the two curves of Fig. 3(b) no entanglement is detected. It is 
surprising that, for the whole region, Qd globally diminishes and tends to be concentrated 
in the im\\-E region, as can be seen in Fig. 3(c). These facts allow one to readily appreciate 
how different is the behavior of entanglement vis-a-vis that of Qd. The role of classical 
correlations can be observed in Fig. 3(d). For the same set of temperatures employed above 
CC decreases i) as we augment T and ii) for increasing values of h, a behavior different from 
that of entanglement: while CC never vanishes, it is larger wherever E = 0. Both E and 
CC coexist for high values of h. We are dealing with a system for which, as we increase 
the temperature, entanglement survives -although barely- for high values of h. This fact 
clearly affects the existence of finite discord- or CC-values. Recall that this was the case 
already at T = 0. The role of the factoring field hf in defining higher or lower values of 
Qd becomes crucial. To further analyze the nontrivial relation between entanglement E 
and quantum discord Qd at finite T it would enlightening to consider a physical system for 
which E would increase with the temperature. Such is the Heisenberg model's scenario, also 
a statistical mechanical model used in the study of critical points and phase transitions of 
magnetic systems [22j E3] . 

III. TWO QUBITS IN THE (FINITE) HEISENBERG MODEL 

First of all note that because of its finitude the system is not immersed in an infinite 
thermal bath. Thus, we cannot stricto-sensu speak of a "temperature" . However, the results 
to be presented are illustrative of the intricacies of entanglement and quantum discord. 
Following the interesting work of Arnesen et al. [23] , we concern ourselves with the issue of 
thermal entanglement but extend the discussion so as to encompass thermal discord. The 
Hamiltonian for the ID Heisenberg spin chain with a magnetic field of intensity B along the 
z-axis reads 



where u l x stand for the Pauli matrices associated to the spin i. Periodic boundary con- 



N 




(14) 



ditions are imposed (cr^ 



JV+l 



aV). Jh is the strength of the spin-spin repulsive interaction 
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(only the anti-ferromagnetic (J# > 0) instance is discussed). If we limit ourselves to the 
case N = 2, we deal with two spinors, i.e., with a two-qubits system. So as to speak of 
"thermal equilibrium" we consider the thermal state 



P(T) 



exp( 



H 



k B T> 



Z(T) 



(15) 



with Z(T) the partition function. Expressing both H and p(T) in the computational basis 
|00), |01), |10), |11) we obtain 



H 
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After defining, for convenience's sake, 
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(16) 



2 J H -2B ) 



^wmy 



exp (— 2w — 2y) 



e wp = exp (— 2w) + exp (Qw); 
e wm = exp (-2w) - exp (6w); 
e wpy = exp (-2w + 2y), 

with w = Jn/ksT and y = B/kBT, we also get 
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In this case the concurrence of p(T) reads 

C = 0; 
C 





C'Wp/2 dwm/2 
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for T > T r 



forT < T r 



l + e~ 2 y + e 2 y + e 8 ™' * " ^ 
Recall that there is no entanglement beyond a certain critical temperature T c = 8 Jh/(^b hi 3) 
as can be seen from the previous C— computation. Also, there is a change in the struc- 



ture of the ground state of hamiltonian (14) when the magnetic field reaches the critical 
value B c = 4Jh- In the limit of zero temperature, the ground state of the system may 
be represented by three different pure states: i) for B < B c (non-degenerate), the ther- 
mal state reduces to the singlet state |\& - )(\& - |, ii) at B = B c (two- fold degenerate) 



||^~)(^~| + ||11)(11|, iii) whereas for B < B c (non-degenerate) we have |11)(11|. The pre- 
vious B— distinctions are crucial in order to understand how the concomitant thermal state 
will respond to T— changes. We expect accompanying behaviors from E and Qd whenever 
the initial state is pure (both quantities coincide in such case). Differences should emerge 
for B > B c whenever we study the unexpected behavior of increasing E versus T as far as 
Qd is concerned. The computation of Qd is in the present case analytic and corresponds 
to sin a = ^ for any (3. The pertinent scenario is the subject of Fig. 4 (let us assume 
J H = 1, so that B c — 4). For the 5-range of values that are smaller than the critical value 
B c , entanglement, discord and CC all diminish as the T increases, as shown in Fig. 4(a). 
This behaviour also occurs at B = B c and is depicted in Fig. 4(b). Notice in both cases the 
sudden death of entanglement, whereas the other quantities "survive" in indefinite fashion. 
Fig. 4(c) plots the same quantities for a magnetic field B > B c . Remarkably enough, in 
this case entanglement as well as the quantum discord augment as T increases. In this case, 
again, E suddenly vanishes while the persistence of the quantum discord and CC. [see Fig. 
4(d)] stresses the fact that for B significantly differing from B c , E is minimalfor all T while 
the quantum discord survives. 

Overall, entanglement and quantum discord display similar behaviours -although with 
clear differences- for a finite quantum system [two spins in the Heisenberg model] but become 
radically different from each other when we consider a system in the thermodynamic limit 
(such as the XY-model). 

IV. CONCLUSIONS 

We have compared entanglement E and quantum discord Qd for magnetic systems at 
finite temperatures, comparing their behavior with that of classical correlations as well. It 
is clear that, some similarities notwithstanding, E and Qd behave in quite different fashion 
in the thermodynamic limit. The distinction we are trying to establish here is blurred in 
the case of finite systems. We conclude that for realistic systems E and Qd should both be 
studied in independent fashion, as they reflect on differen aspects of the quantum world. 
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FIG. 1. Color online) Plot of quantum discord Qd (upper solid curve) and entanglement of for- 
mation E (lower dashed curve) vs the external magnetic field h for two qubits in infinite the XY 
model (nearest neighbors), with anisotropy 7 = ^ at T=0. The region around the factorizing field 
hf concentrates maximum Qd. Inset depicts the corresponding classical correlations CC vs h. See 
text for details. 
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FIG. 2. (Color online) Plot of Qd for the same settings as in Fig. 1 for different relative distances 
R = 1 to 2, 3 and co between spins. The further they are separated, the more they collapse into 
a single curve, which is zero for h > 1. Notice that E rapidly tends to zero for all h in the limit 
R — > oo, while the corresponding Qd remain finite. A similar behavior occurs for CC as depicted 
in the inset. See text for details. 
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FIG. 3. (Color online) (a) Value for E vs h for finite temperatures T=0. 01, 0.1, 0.3, 0.5,1 (from top 
to bottom) for R = 1 and 1 = \- Notice how the region of null entanglement spreads form a 
point (at the factorising field hf = yl — 7 2 ) to a region, (b) Plot of the aforementioned region 
of zero entanglement. The upper and lower curves define de limits of h for a given T where null 
E is found. This figure resembles a phase diagram-like plot where the regions of zero and nonzero 
entanglement are defined, (c) Qd exhibits a particular behavior as T increases which tend to be 
maximum within the limits of zero E. (d) CC vs h plot for the same temperatures. An overall 
decreasing tendency is apparent. See text for details. 
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FIG. 4. (Color online) (a) E (solid line), Qd (long-dashed line) and CC (short-dashed line) vs T 
plots for the thermal state of two qubits in the Heisenberg model at the magnetic field B = 1 
(B < B c ). (b) Identical plot for B = B c = 4. (c) Identical plot for B = 8 (B > B c ). (d) Plot of 
the previous quantities for a high value of B (B = 25). See text for details. 
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